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Abstract 

The scalar heavy quark content of the nucleon, Mq = {N\mqqq\N), is relevant 
for computing the interaction of dark matter candidates with ordinary matter, 
while Ms is important for predicting the properties of dense matter. We 
compute A4q in perturbative QCD to 0{a^). As one goes from Mt to M.c the 
leading order contribution decreases as the number of light quarks is dropping, 
while the radiative corrections grow and are all positive. The leading source 
of uncertainty in the calculation is due to the poorly known value of Ais- A 
related calculation suggests that a large value for Ms may be reasonable. 

I. INTRODUCTION AND RESULTS 

Knowledge of matrix elements Aiq = {N\mqqq\N) provides us with some insight into the 
flavor properties of the nucleon wave function. Since scalars typically couple to rriqqq, values 
of the quark mass operator matrix elements are also relevant for dark matter searches, as 
well as for experiments searching for new weak forces. The value of Adg is important for 
kaon condensation in dense nuclear matter [1], and its value has been the subject of a debate 
[2,3]. 

Heavy flavor A^g's have been calculated in perturbative QCD to the leading order in a 
[4]. Using a technique similar to the ones developed in [4-6] and making use of the four loop 
beta function, anomalous mass dimension, and three loop heavy flavor threshold matching 
coefficients available in the literature [7], in Section 2 we extend the calculation of heavy 
flavor MqS to 0{a^). We flnd 
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Mc^-^il + 0.108 + 0.025 - 0.002)M(1 - x^ds) ^ 78.61(1 - x^ds) MeV , 

= -1(1 + 0.071 + 0.009 - 0.002)M(1 - x^ds) = 80.92(1 - Xuds) MeV , 
lb 

Mt = ^{l + 0.0402 + 0.0009 - 0.0007)M(1 - Xuds) = 84.89(1 - Xuds) MeV , (1) 
where M is the nucleon mass (we use M — 938.3 MeV ) , and 

Xuds ~ ^ ' ^ ) 

with A^M being the shift in the nucleon mass which would result from setting = 0, 

dM 



nis 



A,M = / dy 

onis 







ms=y 







dy {N\ss\N)U=y (3) 



and 



From (1) wc see that there are two competing effects as one goes from the top to the charm 
quark matrix element: since the number of light quarks is dropping, the leading order 
contribution decreases, while the radiative corrections grow and are all positive. We also 
note that the greatest source of uncertainty in the above results is the poor knowledge of 
the strange matrix element entering A^M. 

While the same technique is not applicable to Ais directly, since nis ~ A, the scale 
where QCD becomes strong, in Section 3 we calculate the strange quark matrix element in 
the QCD with a hypothetical heavy strange quark and massless u and d quarks and find 
indications favoring a large value of Als for the physical strange quark mass. 

II. THE SIX FLAVOR CALCULATION 

Wc consider six flavor QCD {uf = 6) with three heavy quarks (c, b, t) and three light 
quarks {u, d, s). The running heavy quark masses are defined to be fhh{lj) in the MS scheme, 
and we define the heavy quark mass ruh to be the "scale invariant mass" , namely the solution 
to the equation mh{mh) — ruh- 

For jj, > ruh the theory is described by the QCD Lagrangian with nj "active" flavors; 
at the scale fi = rrih we integrate out the heavy quark h, and until the next heavy flavor 
threshold the effective theory consists of n/ — 1 flavor QCD, with a shifted gauge coupling, 
and a tower of nonrenormalizable interactions suppressed by powers of m|. The O (1/m^) 
contribution, for example, comes from the operator [D ^jG'^^Y / m\ in the effective theory [4]. 
The value of the coupling and of the coefficients of the nonrenormalizable interactions are 
fixed in perturbation theory by matching S-matrix elements in the full and effective theories. 
The effective theory is asymptotically free, and gets strong at the scale A. 

The mass of the nucleon M is given by an expansion in m^^ and mi of the form 
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FIG. 1. Shown here is a schematic plot of the QCD couphng a(/i)^^ versus In /i, in the MS 
scheme. The threshold corrections have been exaggerated for illustrative purposes. Keeping /xq, uq 
and a\ fixed and varying 'mq{fj,Q) one varies the value of A as well the locations of (some of) the 
heavy flavor thresholds mfi{mfi) = m/j. 



M = A X 



n=l h \^hj m=ln=Oh,l ^ ^ \^h, 



(5) 



where h is a generic index for the heavy quarks, / is a generic index for u and d, and A^M 
is defined in (3). The expansion coefficients depend on the details of the dynamics of the 
theory and may depend logarithmically on the ratio A^/m| and/or mi/ A. We perform an 
expansion in mi /A since chiral perturbation theory works well for u and d quarks. The 
rationale for treating the strange quark differently than the light or heavy flavors is that we 
do not want to assume that either the heavy quark or chiral expansions converge well for 
the strange quark. 

We define a fixed reference scale /iq > m,t with a{fio) = cto- The scale A is defined 
by keeping a (A) = a a fixed; the value itself is unimportant to us, only its logarithmic 
derivative matters. Holding gauge couplings and a\ fixed, we can invoke the Feynman- 
Hellman theorem to compute the qq matrix element in the nucleon by differentiating M 
with respect to mq(/io): 



dM 



d 



"0,"A 



■{N\Hqcd\N) 



{N\mq{fi)qq\N). 



(6) 



The equality above follows from the Feynman-Hellman theorem since enters the QCD 
Hamiltonian in the combination m^gg. We have not specified a renormalization scale in the 
matrix elements as they are renormalization scale invariant. 

Let us first consider the chiral limit = = mg = 0. Approximating (5) with 



M = coA + C 



A3 



mt 



(7) 
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we can express the heavy matrix element as 

■ A3 



j^^^M^hi^o) dA 



A dmq{iio) 



l + O 



Mml 



(8) 



Thus to leading order in A^/Afm^ the heavy quark matrix clement is determined by 
how A varies when we vary the "bare" quark mass rhg{no), keeping a\ and the "bare" 
gauge coupling ckq fixed (see Fig. 1). This calculation requires knowledge of the anomalous 
mass dimension, the beta function for the running of the gauge coupling, and the matching 
functions that relate the gauge coupling and the running quark masses in the full Uf and 
the effective {uf — l)-flavor gauge theories. Here we use the notation and results compiled in 
the paper by Chetyrkin, Kuhn, and Steinhauser [7] to which we refer reader for the formulas 
we use and the references to the original results. 

The coupling constant a(/x) = a^'^f\iJ,)/n in the theory with Uf "active" flavors runs 
with the scale /i according to the equation 

(9) 

where the QCD beta function, 

fe>0 

has been calculated to four loop order and the coefficients (3^^^^ are known for A; = 0, 1, 2, 3. 
The current quark mass mq''^\ii) of the theory with n/ active fiavors obeys 



1^' 



d 



mW)(/x)=m("/)(/x)7^)(a), (10) 



where the anomalous mass dimension, 

fc>0 

has been calculated to four loop order and the coefficients 7^1^ ^-^e known for A; = 0, 1, 2, 3. 
At a heavy quark threshold n = ruh the couplings and the quark running masses in the 
"full" and the effective theories are related by the matching conditions 

a("/-^)(^/.) = a("/)(m,)Cj(a("/)(^^)), 

m^q'-'\mh) = m^^^\mh) Ua^''^\mh)). (11) 

The functions (^g and Cm have been computed in perturbation theory to three loop order. 

Here a^"f\mh),fn"q^^\mh) and a^"'f~^\mh),mq'^ ^\fTih) are the coupling and quark masses 
evaluated just above and below the flavor threshold, respectively. 
Equations (9) and (10) may be integrated to give 
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ln[ 



rrih 



da— 7 — T——. 

a{mn) (a) 



(13) 



In order to obtain the desired quantity in (8) we differentiate (12) witli respect to 
In rnqdio), keeping /iq, ao, and cka fixed and taking into account tfie flavor tfiresfiold discon- 
tinuities. Tfie result is 



D 111 A 
din fhgiiio) 



El- 



pi-^f ^){Q{anfh)anfh) da^fh J d In mgino) 



(14) 



Here Un^h = a^^f^rrih) and the sum is over the heavy flavors lighter than the flavor q and 
the q itself. To calculate 91nm^/91n rhq{fiQ) we differentiate (13) with respect to In mg(/xo) 
keeping /iq, ckq fixed and taking into account the flavor threshold discontinuities in the 
running quark masses and coupling constant. 

Now let us take into account the light quark masses. Working in the isospin limit 

= ''TT'd = we approximate (5) with 



M = coA + 2 cio/,i + A,M + O 



A3 



mi 



m 
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and the heavy matrix element expression is modified to 



Mq = M{l-Xuds) 



A drugdio) 



l + O 



A 



A3 



mt, 



Mml' MA^ 



(15) 



(16) 



where x^ds^^ defined in (2); E^riv — 45 MeVis defined in (4) and is known from pion-nucleon 
scattering, while the value of the strange matrix element, Ms — 200 MeV , that enters A^M 
from (3) is quite uncertain and varies by hundreds of MeV depending on the method of 
calculation [2,3]. 

Wc used the following values of parameters: Mz = 91.18 GeV , a{Mz) = 0.117, rric = 1.2 
GcV, nib = 4.2 GeV and mt = 175 GeV [8]. The numerical results obtained are shown in 

(!)• 



III. CALCULATION IN THE THREE FLAVOR TOY MODEL 

The value of Ais has been calculated in baryon chiral perturbation theory. The leading 
order calculation predicts linear growth of the matrix element as a function of m^, however, 
at one loop order one gets sizable negative contributions from the meson loops proportional 
to m^/^ [2,3]. This indicates that chiral perturbation theory may not be the right tool for 
the calculation. Calculations of Ms in the Skyrme model lead to the same conclusion [9]. 
This is unfortunate since the value of Ms may have important consequences for the nature 
of dense matter [1], and the fate of supernova remnants [10]. It could also be relevant for 
the properties of the strongly bound kaonic system K^ppn observed recently [11]. 

One may, however, gain qualitative insight into this problem using a rather different 
approach. We calculate the strange matrix element in three flavor QCD with a hypothetical 
heavy strange quark (such that ^ A) and massless u and d quarks using the technique 
of the previous section. 
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FIG. 2. Shown here is the perturbative result for A4s/M = {N\msSs\N) /M versus in 
GeV calculated in the toy three flavor model with massless u and d quarks and a heavy s quark. 
At the leading order Mg/M = 2/29 ~ 0.069 (the red line); the green line corresponds to the NLO 
result, the blue line - NNLO and the black dashed line - NNNLO. Perturbation theory is clearly 
breaking down for rUg < 1 GeV . 



In the effective nj = 2 theory the mass of the nucleoli M is given by an expansion in 
m7^ of the form 



M 



n=0 



(17) 



where the coefficients c„ will depend logarithmically on the ratio A/m^. Defining a reference 
scale fiQ > nis with fhs{fio) = rriso and gauge coupling a(//o) = Q^o, and using the same line 
of argument as in the previous section, we can express the strange matrix element as 



rriso dA 



A dm 



sO 



l + O 



A3 



(18) 



Thus to leading order in A^/Mm^ the strange matrix element is determined by how A varies 
when we vary the "bare" strange quark mass, keeping and the "bare" gauge coupling ao 
fixed. 

In order to obtain the desired quantity in (18) we differentiate both sides of (12) and 
(13) with respect to Inmgo, keeping fiQ, ao, and fixed and taking into account the strange 
quark threshold discontinuity. We get 



{N\msSs\N) 



(9 In A 
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27^ V) 



/3(3)(( 



d 



/?(2)(C2(a) a) da 
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FIG. 3. 



where a = a^^\ms). The couphng a^^^rris) is calculated by running down from /j, = Mz to 
li — nis using (12) with the real world initial condition a{Mz) = 0.117. 

Fig. 2 shows Ms/M = {N\msSs\N) / M as a function of nis calculated through four loop 
order in perturbation theory. It indicates that QCD may favor a bigger value of the strange 
matrix element for the physical m^. A reasonable guess for Ms{fns) in the nonperturbative 
regime consistent with what we know about Ms for the values of rUs, both larger and smaller 
than its physical value, is sketched in Fig. 3. 
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